Can the local energy— momentum conservation laws be derived solely from field 
equations? 
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Abstract. The vanishing of the divergence of the matter stress— energy tensor for General Relativity is a particular case 
of a general identity, which follows from the covariance of the matter Lagrangian in much the same way as (generalized) 
Bianchi identities follow from the covariance of the purely gravitational Lagrangian. This identity, holding for any covariant 
theory of gravitating matter, relates the divergence of the stress tensor with a combination of the field equations and their 
derivatives. One could thus wonder if, according to a recent suggestion [1], the energy-momentum tensor for gravitating 
fields can be computed through a suitable rearrangement of the matter field equations, without relying on the variational 
definition. We show that this can be done only in particular cases, while in general it leads to ambiguities and possibly to 
wrong results. Moreover, in nontrivial cases the computations turn out to be more difficult than the standard variational 
technique. 

In a recent paper [1] Accioly et al. have observed that for some weU known cases of classical fields 
interacting with gravitation (e.g. scalar field, electromagnetic field), upon contracting the dynamical equation 
for the matter field with a suitable linear combination of covariant derivatives of this field, the resulting 
expression represents the vanishing of the covariant 4-divergence of some rank-two tensor. In the cases 
considered in [1] , the latter turns out to coincide exactly with the stress-energy tensor of the matter field, 
according to the usual variational definition. 

Let us explain the geometrical origin of this phenomenon. For the reader's convenience, we recall the 
derivation of the strong conservation law (sometimes called "Bianchi identity for matter" ) in the case of metric 
theories of gravity, in the form suitable for our aim. We denote arbitrary matter fields by V'a; the signature 
of the metric 5^1^ is chosen to be ( — h ++); the expressions V^i/ and both denote covariant derivation 
relative to the metric g^i, and we set c = SttG = 1. In accordance with [1], we assume that pure gravitation 
is described by the usual Einstein-Hilbert Lagrangian, although the considerations below hold, mutatis 
mutandis, under more general assumptions. The matter Lagrangian is L{g, tp) = L{g^,j, i?^^,^., i/;^, 'ipA■,^l)'■ the 
dependence on g includes a possible non-minimal coupling to the curvature. The resulting action is 



The (Hilbert) stress-energy tensor and the l.h.s. of the equation of motion for the field V'A are defined, 
respectively, as 



The action (1) is coordinate invariant, so it is left unchanged by any infinitesimal point transformation 
^ x^^ -\- {x) , with being an infinitesimal vector field [2-4] . For simplicity, we assume that = 
on the boundary of the integration domain fJ, so that the latter is mapped onto itself by the infinitesimal 
transformation. The variation of the metric is 
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and 
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where £5 denotes the Lie derivative. Similarly, for the matter variables the variation Si/ja = — -Cj^a can be 
expressed [4] in terms of the covariant derivatives of 'ipA, which have the general form 



M/3 
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where the coefficients Z^^a linear functions of ip and depend on its tensorial rank^; namely, one has 
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The invariance of the action implies 



= bS=\ ^ 

Jn 



(6) 



plus a surface integral which vanishes under our assumptions on Using (3) and (5), dropping again a 
total divergence and taking into accoimt the (metric) Bianchi identity VG/^jy = 0, the integrand in (6) 



becomes ^^[V'-'T^^ - E'^V^iPa - Vf){E-^ZA 



The vector field being arbitrary in the interior of f2, the 



vanishing of the integral (6) entails the following identity: 



(7) 



which implies the local conservation law VT^u = for solutions of the field equation E^{il),g) = 0. 

Equation (7), which is quoted in this form in [4] but was essentially contained in earlier work (see 
e.g. Ref. [2]), shows that contracting the free index of the l.h.s. E^ of the matter field equation with V^tA/I 
always yields the full divergence of a rank-two tensor. This is a universal property, holding for any matter 
and any (generally covariant) Lagrangian. 

Apparently, (7) provides a possible way to compute the tensor T^^ directly from the field equations, 
and therefore to obtain the r.h.s. of the Einstein equation without having to deal with the action principle. 
This possibility has been advocated by the authors of [1], who however identify directly the stress-energy 
tensor with the full expression occurring in the divergence on the r.h.s. of (7), thus overlooking the term 
S/^^E'^ZAvfi)- Of course, this term vanishes for solutions of the field equation E^ = 0, but for most purposes 
the correct stress-energy tensor of the theory should be unambiguously defined for any field configuration, 
and not only for exact solutions. Only for the examples considered in [1] can the correct stress-energy tensor 
be obtained according to their prescription - in the first two examples (scalar massive Klein-Gordon field 
and non-minimally coupled scalar field), just because for a scalar field this term vanishes. For the subsequent 
example, i.e. for the electromagnetic potential coupled (either minimally or non-minimally) to gravity, 
the identity (7) becomes 

E"A,,^ = y"{T^, + E,A^). (8) 
If one multiplies Maxwell equations F'^'^ = by A^-a, after a number of rearrangements one arrives at 



R^'^^-'F^o^A^, where r 



■I"" = -Ac,'''F''" + ^g'"'Fo,f)F°''^ 



is the canonical energy-momentum tensor, which is non-symmetric and gauge-dependent and as such has 
no physical meaning. This may bo particularly confusing if one makes this computation in flat space, where 
this tensor is conserved. However, the authors of [1] multiply instead the electromagnetic equation E" = 
by Ff^,y = VfiAi, — VvAfi, thus producing an additional term which cancels exactly the last term on the 
r.h.s. of (8): in fact, E"Ff^„ = E" A,,.,^ - V,y(£"'A^) + A^y „E" and the divergence of the field equation, 
Vi/£''', vanishes identically in the case discussed. Unfortunately, in general it is not possible to get rid of the 
additional term in (7) by contracting the field equation with a suitably modified combination of covariant 



^ For instance, if (j}A is a collection of scalar fields then Za^ ce = 0) if deals with a vector field V''' one 
has Z^^^a = V'"'^^) and so on. 



derivatives of the field: namely, the prescription fails to work whenever V„E" does not vanish identically. 
For example, consider a vector field with a Lagrangian 

X = S'^^S'^^ where S^^ = A^.^ + A^.^ ^ S^^. (9) 

The Eulcr-Lagrangc field equations are Maxwell-Iike, S^'^ -j^ = 0, but now S'^'^-f^,^ docs not vanish for arbitrary 
A^. Thus, the r.h.s. of (8) includes a term, A^S^'^-x^, which contains hnearly the third derivatives of the field 
Afj, and therefore cannot be cancelled by any linear combination of the form c^^Aa-fjE^ (with c"^ constant 
coefficients) which would contain only the first and second derivatives. 

We further remark that any manipulation of the equations along the ideas presented in [1], aimed at 
producing the correct result in at least a reasonably wide class of theories, is likely to exploit the freedom of 
multiplying the field equations by a constant: then, the resulting stress-energy tensor would be determined 
up to a multiplicative factor (of any sign). However, to identify T^^ in the Einstein field equation with 
the physical energy and momentum density, the correct numerical factors (which are suppressed in many 
papers) should be included in matter Lagrangians. This is the most elementary reason why the information 
contained in the Lagrangian cannot be fully replaced by the knowledge of the matter field equations. The 
actual risk inherent in playing with the field equations alone is well illustrated by a computation given in [1] 
for the no n- minimally coupled scalar field. The minimal coupling, according to their equation (8)[1], should 
be recovered when /(</>) reduces to any constant. However, the terms due to the non-minimal coupling 
in the stress tensor given on p. 1165 disappear only if /((^) = 0. In this case, however, the stress-energy 
tensor reduces to exactly half of the expression (7)[1] given for the minimally coupled field. Setting instead 
/(0) = 1/(2k), as is suggested in [1], leads to a wrong result: the Einstein equation becomes the correct one 
(multiplied by ^), but the expression for T^^" equals (7)[1] only for solutions of the Einstein equation. 

We shall discuss below whether a more refined prescription might be devised, taking into account the 
additional term E^Zauh in (7), to implement effectively the proposal of [1]; before that, let us address a 
side aspect which may also cause misunderstanding. 

The authors of [1] make a preliminary distinction between the theories involving only minimal coupling 
and those including non-minimal coupling. In doing so, they assume that minimal coupling occurs if the 
matter Lagrangian contains only the metric g and not the Levi-Civita connection F, which is incorrect. The 
minimal coupling rule states that in a freely falling local reference frame (i.e. in a locally geodesic coordinate 
system) the matter Lagrangian should reduce to the form of the flat-space matter Lagrangian. This implies 
that a minimally coupled matter Lagrangian should not contain curvature terms; however, it will necessarily 
contain first derivatives of the matter fields, and in general this entails the occurrence of Christoffel symbols 
to ensure covariance. In the cases of scalar fields and of electromagnetism it is possible to get rid of the 
Christoffel connection in the matter Lagrangian, but in other cases a "minimally coupled Lagrangian" not 
including F does not exist at all. Examples arc provided by spin-^ Dirac field [5], linear spin-two fields in a 
Ricci-flat background [3], spin-| field [6,7] and the vector field (9) considered above. It is worth mentioning 
here that there might exist cases where the metric connection F occurs in the matter Lagrangian not only 
through covariant derivatives of the fields. Consider for instance a matter field 7"^ which transforms as 
a (non-metric) linear connection. Then one defines a tensor Q'^^ = 7^^, — V^^ and expresses the matter 
Lagrangians in terms of Q and its covariant derivatives (with respect to F). Should such a coupling be 
considered a non minimal one? 

In any case, we stress that whether the matter fields are coupled to gravity minimally or non-minimally 
is irrelevant by itself in this context. What really affects the discussion is whether the full Lagrangian can be 
split into a pure gravitational part plus a matter term, or not; if it cannot, the stress energy tensor itself is 
not well-defined. Now, non-minimal coupling does not necessarily prevent such a splitting of the Lagrangian. 
However, for some non-minimally coupled theories one is led to define an effective stress-energy tensor (for 
instance in scalar tensor tcories; see Ref. [8], eq. (2.25) (2.28)^) which does not correspond to a separate 
matter Lagrangian; as long as one looks at the gravitational equations, this tensor is indistinguishable from 
a genuine variational stress-energy tensor, yet in this case trying to derive the proper conservation laws from 
the field equations would lead to ambiguities. 



^ In [8], eq. (2.25) is affected by a misprint and should read G/jtv = ■■ ■ = 9^^{(t>,g) + ^T^^{^,g). 



A different, but not unrelated, problem is connected with the fact that while dealing with gravitational 
theories including non-minimal coupling, possible redefinitions of the fields (including the metric) are often 
considered in the literature (see for instance [9], where a particular application of the general method 
described in [10] is discussed in detail); this causes further ambiguities, on completely different grounds, 
on the definition of the correct physical stress-energy tensor [11-13]. Such ambiguities can be removed by 
a careful analysis of the Lagrangian formulation of the model [8], while a prescription based on the field 
equations alone would be unsuitable for a rigorous approach to this problem. 

Let us finally come back to the main question: docs the universal formula (7) allow us to obtain a 
practical method to compute the stress-energy tensor? The identity (7) allows to single out T^j/ only up 
to the addition of an arbitrary tensor being identically divergenceless, for instance V^"^ = VaW^""', where 
y^n^i^a _ yyltii^a] jg g^j^y totally antisymmctric tensor (then Ricci and Bianchi identities imply 'Vi^V'^'^ = 0). 
The procedure should then at least be supplemented by the requirement that the resulting rank-two tensor 
be symmetric; yet, in general it is unclear whether this would ensure uniqueness. The uniqueness cannot 
be restored by requiring that the kinetic part of T^^ be quadratic in the field variables, because for a non- 
minimal coupling in general the variational stress tensor contains terms linear in the highest derivatives [8] . 
Furthermore, the method remains highly non-algorithmic, not because the manipulation to be performed on 
the field equation would be only vaguely defined (as it appears in Ref. [1]) but rather because recasting the 
l.h.s of (7) into the divergence of a symmetric tensor needs several non-trivial tricks to be found ad hoc. 

For instance, the computation for the vector field occurring in (9) using (8) would be quite tricky; it is 
much easier to use a generic formula for A/j^ following from (2), 



where Qf"" = dL/dA^.^. 

For a spin-two field, which is represented by a symmetric rank-two tensor the identity (7) reads 



and E^^"^ is so complicated [3] that guessing the appropriate manipulations to get from (11) would hardly 
be successful. 

We conclude that to compute the stress-energy tensor it is in general both safer and easier to rely on 
the variational definition (2), which provides an algorithmic and unique prescription. The computation can 
be performed almost straightforwardly using standard computer packages for tensor calculus (see e.g. [14]). 
On the other hand, if the procedure suggested in [1], rather than a purely computational trick, is intended 
to provide an alternative way to define the matter source term in a general-relativistic gravitational theory, 
circumventing the need to introduce a Lagrangian density for the model, then the whole approach is 
misleading. In fact, the definition itself of stress-energy tensor, as well as the property that the matter 
field equation can be recast into a full divergence by a suitable manipulation, rely on the existence of an 
action principle from which both the gravitational (e.g. Einstein's) equation and the matter field equation 
should be derived. 
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T^.{A,g) = Lg^, - 2^ + 2V^ [q^^^^ A^^g^^^) - Va , 



(10) 




(11) 
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